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1. INTRODUCTION 

Due to the severe environment problems resulted from using fossil fuels, it is very important to 
propose clean energy techniques in order to keep human sustainability. Previous fact gives the modern focus 
on nuclear energy generation [1], [2]. It was consedered that fission energy is a pivotal clean energy that does 
not emit greenhouse gases and could replace fossil fuels in an efficient way with huge amount and economic 
competitiveness [3], [4]. Because the pressurized water reactor (PWR) is the most commonly applied nuclear 
reactor, safe and stable operation of the PWRs is considered very important issue for the nuclear industry [5], 
[6]. Power-level regulation is important in order to maintain a satisfactory operation of PWRs. The main rule 
of power level control is to determine the injection and withdrawal speed for the reactor controlling rods to 
maintain the generated power at a required rate [7], [8]. 

The most widespread and simple control element in different fields of industry is the proportional- 
integral-derivative (PID) controller because of its singularity in its plain structure and easy implementation in 
modern applications and complex industries [9], [10]. Although the PID controller is very simple in 
construction and implementation, it needs uninterrupted refinement in order to enhance its performance [11]— 
[13]. The PID performance can be further enhanced through the process of making use of fractional order 
derivative and integral. The two added parameters in the FOPID controllers give more opportunities for 
flexibility, control and stability [14]. 

In general, the tuning methods for FOPID controller can be classified into analytical, numerical and 
rule-based ones [15], [16]. The controller parameters can be analytically determined by solving nonlinear 
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equations in order to achieve the required specifications in the gain and phase margin and crossover 
frequency. But this method is applicable only when the equations are small in number and plain. So, it is 
difficult to get FOPID parameters for the the PWR. For the rule-based method, it can determine the required 
parameters depending on empirical tuning rules. However, the consedered system usually must be a system 
that having an S-shaped step response [17]. 

On the other hand, the numerical method - which is usually an optimization-based method - is a 
relative suitable solution for tuning FOPID for the PWR [18], [19]. Three optimization techniques is used in 
this study which are particle swarm optimization (PSO) [20], [21], gray wolf optimization (GWO) [22], [23] 
and ant lion optimization (ALO) [24], [25] and the evaluation will be based on the integral square error (ISE). 


2. PROBLEM FORMULATION 

The nominal PWR model used in FOPID controller designing is from order five. It is point kinetics 
having only one delayed neutron group. It also has a temperature feed-back from agglomerated fuel and 
coolant temperature computation that can be listed as follows [6], [8], [26], [27]: 


Mr — eE n, +Ec, (1) 
a = An; — Ac; (2) 
Pease peg it ate (3) 
an z a kg a1, mun T, + au T, i 
acer = G2, (5) 
dp = 6p, + ar (Ty To) + ecto) es seeto) (6) 


where: 

Poa =Rated power level (MW) 

n, = Neutron density relative to density at rated condition 

fy = Portion of reactor power deposied in fuel 

Hy = Fuel heat capacity 

Hc = Coolant heat capacity 

Q = Fuel and coolant heat transfer coefficient 

M = Mass flow rate times coolant heat capacity 

T; = Average temperature of reactor fuel 

T, = Temperture of the coolant leaving the reactor 

Te = Temperture of the coolant entering the reactor 

a, = Fuel temperature reactivity coefficient 

a, = Coolant temperature reactivity coefficient 

Tyo = Equlibrium average reactor fuel temperature 

Tio = Equlibrium temperature of the coolant leaving the reactor 
Teo = Equlibrium temperature of the coolant entering the reactor 
ôp = Reactivity 

6p, = Reactivity because of the control rod 

Z, = Control input, control rod speed (fraction of core) length per second 
G, = Total reactivity worth of controlling rod 


The PWR used in this study can be modeled as 5“ order nonlinear model having only delayed neutron group 
and pair of thermal feedback system. Table | shows the parameters of the PWR. 
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Table 1. The parameters of the PWR [28] 


Parameter Value Unit 
Poa 2500 MW 
ff 0.92 
Ug 26.3 MW. s/°C 
Ue 71.8 MW. s/°C 
Q 6.6 MW/°C 
M 102 MW/°C 


3. METHOD 

In this section, the three optimization techniques are used to tune the five parameters of the FOPID 
controller applied to the process (PWR model) and the selected objective function to be minimized is the 
integral square error (ISE) for a unit step response. The simulation is performed using MATLAB software. 
Figure 1 shows the generalized block diagram for the FOPID controller with the process. The transfer 
function of the FOPID controller has the following form: 


G(s) = K, + í + Kas” (7) 
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Figure1. Block diagram for the FOPID controller with the process (PWR model) 


The PSO technique relys on the intelligent motion rule of a swarm that proceeds within a search 
area, searching for the best solution to it [29]. So, every particle within vacant region is treated as an N- 
dimensional point that trys to optimize its posture with the use of the current location and the speed [21]. The 
GWO is a modern inspired algorithm that simulates the directing pyramid and hunting method of grey 
wolves in real life. 

This hunting process of grey wolves consists of three steps: looking for victim, encircling the 
victim, and attacking the victim. Grey wolves normally move in groups and are divided into four groups. The 
first group wolves are the commander, the second group wolves support the first group in obligations and can 
replace them if they die, the third group are the hunters, keepers and explorers of the group, the forth group 
refers to the young group members [22]. The ALO algorithm simulates the hunting mechanism of antlions in 
real life. In ALO, a selected random walk is usedin order to stochastically represent the ant’s motion toward 
the nourishment place. During optimization procedure, ants modify their locations in each iteration within 
walk sat stochastically determined. every random walk step is calculated using the min-max normalization 
technique to confirm the rate of random walks within the trait space [25], [30]. 

Table 2 shows the values of the parameters of the FOPID controller and the integral square error 
when the three optimization techniques are used to minimize the ISE due to unit step variation in the 
reference power. It is clear that the ant lion optimization technique gives the best ISE. 


Table 2. The parameters of FOPID controller and the ISE due to unit step response 
Variable Ky K; Ka À H ISE 
PSO 70 25 17.565 0.679 0.4995 0.002331 
GWO 68.23 6.3748 7.2609 0.7048 0.0138 0.11017 
ALO 69.9876 24.9931 17.9227 0.6345 0.5 0.002328 
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4. TESTS AND SIMULATION RESULTS 

Several tests are implemented using MATLAB to evaluate the efficiency of the three optimization 
techniques as shown in Figure 2. Figure 3 shows the normalized output power responses due to step variation 
of reference power from 85% to 90% at time equals 25 seconds. Figure 4 shows the normalized output power 
responses due to step variation of reference power from 100% to 90% at time equals 50 seconds. Figure 5 
shows the normalized output power responses due to step variation of reference power from 80% to 85% at 
time equals 20 seconds and then back from 85% to 80% at time equals 40 seconds. In order to qualitatively 
judge the performance of the three optimization techniques on tuning FOPID controller parameters, Table 3 
shows the ISE results and comparison. 
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Figure 2. MATLAB simulink implementation of the FOPID controller for the PWR model 


Normalized Output Power 


25 25.5 26 26.5 27 27.5 28 
Time 


Figure 3. Perunit output power responses due to step variation of reference power from 85% to 90% 
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Figure 4. Perunit output power responses due to step variation of reference power from 100% to 90% 
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Figure 5. Perunit output power responses due to step variation of reference power from 80% to 85% and then 
back from 85% to 80% 


Table 3. The integral square error for the responses in Figure 3-5 


Technique ISE (Figure 3) ISE (Figure 4) ISE (Figure 5) 
FOPID-PSO 0.0020514175 0.002391215 0.001911005 
FOPID-GWO 0.007242511 0.008494478 0.006430943 
FOPID-ALO 0.002045736 0.002213185 0.0018 16844 


It is noted that the PSO based FOPID controller has performance close to the ALO based FOPID 
controller. However, the ALO based controller has lower integral square error. Also, the GWO based FOPID 
controller has the least overshoot but has relatively long settling time and the worst ISE. 


5. CONCLUSION 

Three optimization techniques are used to design a fractional order PID controller to control the 
output generated power of a pressurized water reactor. The simulation results show that the ant ALO 
technique gives the best performance in tracking the reference power based on the integral square error 


index. Also, ALO based FOPID controller gives the minimum settling time with zero steady state error. On 
the other hand, the GWO based FOPID controller has the higher settling time and the worst ISE. 
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